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Abstract 

The  paper  deals  with  a  concept,  mathematical 
foundation  and  results  of  testing  a  prototype  of 
onboard  system  providing  an  aid  for  pilot's 
control  actions  during  more  or  less  long-term 
manoeuvres,  such  as  takeoff/climbing,  waypoint 
flight  through  the  threat  zones,  surface-based 
and  air  target  attack,  descent  and  landing  ap¬ 
proach.  The  key  issue  of  this  system  is  to  pro¬ 
vide  an  algorithm  allowing  on-line  prototyping 
of  spatial  trajectories  for  these  manoeuvres  and 
visualisation  of  the  reference  trajectories  to  a 
pilot  on  3-D  flightpath  display  for  further 
tracking.  As  an  efficient  tool  for  onboard  opti¬ 
misation  the,  present  paper  considers  different 
modifications  of  a  direct  method  of  calculus  of 
variations,  implemented  to  the  tasks  of  flight 
dynamics  first  time  by  Prof.  Taranenko  in  the 
early  60’ s.  The  paper  contains  a  brief  survey  of 
the  methods  used  so  far  by  the  other  researches 
in  the  field,  discloses  main  ideas  of  proposed 
approach,  based  on  independent  optimisation  of 
a  trajectory  and  speed  profile,  shows  how  to 
assure  a  required  convergence  robustness.  The 
paper  is  concluded  with  the  reference  to  the 
semi-natural  and  flight  test  of  the  proposed 
system  prototype,  carried  out  in  Russia  in  mid 
90 ’s  in  Zhukovskiy  Air  Force  Engineering 
Academy  and  Gromov  Flight-Test  Institute. 

1  Introduction 

It  is  well  known  that  a  biocybernetic  cabin  is 
considered  by  aviation  specialists  of  leading 
countries  as  one  of  the  most  progressive  ways  of 
modern  and  perspective  aircraft  (a/c)  efficiency 
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increase  (see,  for  example,  Refs.  [1,2]).  A  hard¬ 
ware  and  software  of  the  last  one  (also  known  as 
a  Pilot’s  Associate  (PA)  or  Electronic  Co-pilot) 
should  provide  a  pilot  with  the  relevant  infor¬ 
mation  and  intelligent  support  during  flight  mis¬ 
sion  fulfilment  -  for  all  tasks  from  flight  plan¬ 
ning  and  preparing  up  to  the  landing  and  post- 
flight  analysis.  Pilots  expect  such  support  would 
increase  the  total  “board  intellect”,  i.e.  summary 
possibilities  of  man-machine  (MM)  system. 
Therefore,  this  system  would  allow  them  to 
solve  the  different  complex  flight-mission’s 
tasks  creatively,  not  as  a  robot.  Furthermore,  15- 
year  a/c  accidents  statistics  [3]  and  pilot’s  inter¬ 
viewing  [4]  strongly  shows  that  because  of  the 
fact  that  specifically  control  errors  evoke  a  ma¬ 
jority  of  the  flight  accidents,  practically  all  pilots 
have  the  urge  requirement  of  their  control  ac¬ 
tions  support.  (By  the  way,  they  consider  this 
last  task  as  the  most  difficult  to  be  supported 
[4].)  The  recent  progress  in  a/c  complexity  and 
potential  possibilities  makes  this  requirement 
highly  urgent.  Thus,  the  question  ‘whether  we 
have  any  means  to  provide  a  pilot  with  a  rele¬ 
vant  cognitive  prompt  to  assist  him  in  control¬ 
ling  the  a/c?’  is  vital.  As  a  matter  of  fact,  this 
question  can  be  split  into  two  separate  questions 
-  the  question  of  cognitive  representation  of 
guidance/control  prompt  and  the  question  of 
relevant  calculation  of  the  reference  controls. 

1.1  Perspective  flightpath  displays 

According  to  pilots  [4]  the  most  desirable 
prompt  during  manoeuvring  would  be  the  refer¬ 
ence  trajectory  (RT)  shown  on  a  head-up  (HUD) 
and/or  head-down  displays  in  a  view  of  a  “road- 
in-the-sky”  (RiS)  (also  being  referred  to  as  a 
tunnel-in-the-sky,  pathway-in-the-sky,  highway- 
in-the-sky,  etc.)  or  in  the  other  words  on  a  3D 
flightpath  display  (FD). 
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Probably  the  first  time,  FD  was  referred  in 
the  technical  literature  in  the  early  50’ s  in  con¬ 
junction  with  the  US  Army  Navy  Instrumenta¬ 
tion  Program  [5].  The  goal  of  this  program  was 
to  develop  an  improved  MM  interface.  Of 
course  at  that  time  because  of  the  lack  of  appro¬ 
priate  computer  technologies  not  many  of  FD’s 
formats  proposed  in  the  context  of  the  men¬ 
tioned  program  and  latter  researches  (by  Hoo¬ 
ver,  Farrand,  Wilckens,  Jensen,  and  others)  was 
implemented  and  flown.  The  major  area  of  in¬ 
terest  at  that  time  was  to  employ  perspective 
FDs  for  the  instrument  landing  guidance. 

Latter  research  (done  by  Wattler,  Mulley, 
Grunwald,  Adams,  Logan,  Jauer,  Quinn,  Mitori, 
et  al.  in  the  80’ s)  addressed  many  issues  of  FD, 
including  RiS  representation:  its  cross-section 
shape  (straight  and  curved  U,  brackets,  tiles, 
square,  rectangle,  etc.),  cross-section  connec¬ 
tions,  sizes  and  so  on.  They  also  discussed  a  RiS 
as  an  alternative  to  the  flight  director  and  pro¬ 
posed  to  employ  a  synthetic  terrain  to  improve  a 
perception  FD’s  format  [6], 

One  of  the  first  FD  prototypes  was  flown 
by  Hoover  in  the  early  80’ s.  A  total  of  90  hours 
of  ground  simulation  and  20  hours  of  in-flight 
operation  demonstrated  that  such  display  pro¬ 
vides  the  pilot  with  adequate  information  to  per¬ 
form  all  types  of  steady  state  manoeuvres  (in¬ 
cluding  takeoff/climbing,  cruise  flight,  approach 
and  landing)  without  reference  to  conventional 
parametric  displays  or  the  real  world.  Based  on 
the  results  of  the  flight  it  was  stated  that  the 
concept  of  continuous  command  information  is 
one  of  the  most  significant  innovations  that  FD 
can  provide  to  MM  interface.  This  interface 
serves  as  external  memory,  eliminating  the  ne¬ 
cessity  for  memorising  each  segment  of  the 
flight  mission  plan  and  referring  to  a  navigation 
chart,  reducing  the  load  on  working  memory  [6]. 
Hoover  et  al.  [7]  indicated  the  operational  ad¬ 
vantages  by  stating  that  no  coaching  by  a  second 
pilot  with  regard  to  each  upcoming  event  was 
necessary  when  the  evaluation  pilot  was  flying 
the  RiS.  With  respect  to  navigational  awareness, 
they  reported  that  while  flying  the  FD,  the  ma¬ 
jority  of  the  pilots  stated  that  'there  was  never 


any  question  as  to  where  they  were,  what  they 
were  doing,  or  what  they  should  be  doing'  that  is 
not  a  case  with  a  common  symbolic  format.  Be¬ 
cause  of  mentioned  promises  FD  was  included 
into  announced  PA  program  as  an  output  inter¬ 
face  in  the  mid  80’ s  (Fig.l)  [1], 


In  the  90’ s,  the  research  teams  led  by 
Grunwald,  Swenson,  Theunissen,  and  Sachs 
were  able  not  only  to  proceed  with  semi-natural 
experiments  using  the  modern  powerful  com¬ 
puters,  but  also  to  flight  test  their  FDs  on  the 
different-type  a/c.  They  reported  a  significant 
improvement  in  pilot  situational  awareness  and 
mission  effectiveness,  as  well  as  a  decrease  in 
training  and  proficiency  time  required  for  a  near 
terrain,  night-time,  adverse  weather  flight.  They 
also  studied  an  implementation  of  additional 
‘dimension’  -  superimposed  predicted  vehicle 
position  (4D  approach)  [6].  Fig.2  shows  pilots’ 
desire  to  have  a  RiS  prompt  at  the  different 
stages  of  flight  [4], 


0%  25%  50%  75%  100% 


Fig.2  Pilots’  desire  to  use  FD  at  different  stages  of 
flight 


All  preceding  research  were  based  on  some 
pre-set  trajectory  calculated  beforehand  off-line 
(usually  it  was  a  straight  glide-path  trajectory  or 
piecewise  trajectory  consisting  of  straight  lines 
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and  pieces  of  spirals  or  straight  terrain¬ 
following  trajectory).  Today’s  almost  1GHz- 
frequency  computers  with  enormous  storage  ca¬ 
pacities  and  advanced  software  allow  designing 
any  complex  adaptive  FD  and  its  usage  on 
board.  However,  one  obstacle  still  prevents  the 
wide  employment  of  these  displays  -  a  lack  of 
reliable  method  allowing  on-line  prototyping  of 
any  spatial  trajectory  for  any  stage  of  flight. 

1.2  Pilot’s  control  actions  support 

As  mentioned  above,  we  have  to  be  able  to  cal¬ 
culate  optimal  trajectories  for  a  majority  of  ma¬ 
noeuvres  in  accelerated  scale  of  time.  That  is  a 
main  goal  of  the  proposed  on-board  pilot's  con¬ 
trol  actions  support  system  (PSS),  considered  as 
an  element  of  supervising  PA  (Fig.3)  [8].  The 
structure  of  such  PA’s  subsystem  for  the  major¬ 
ity  of  different  tasks  to  be  supported  is  shown  on 
Fig.4.  In  real  flight  for  the  real  strategic  and/or 
tactic  conditions  for  real  state  of  the  a/c,  its  sub¬ 
systems  and  even  the  pilot,  PSS  provides  com¬ 
putation  of  RT  to  be  visualised  on  FD. 

Since  the  requirement  of  optimisation  ra¬ 
pidity  and  analytical  representation  of  the  ex¬ 
tremal  solution  are  the  main,  an  employment  of 
classical  methods  of  calculus  of  variation  be¬ 
comes  strongly  limited.  This  has  led  to  a  search 
for  variational  methods  of  a  different  kind, 
known  as  direct  methods  (DMs)  (which  do  not 
entail  the  reduction  of  variational  problems  to 
problems  involving  differential  equations).  The 
fundamental  idea  of  DMs  is  to  consider  a  varia¬ 
tional  problem  as  a  limit  problem  of  the  extreme 
of  a  function  of  a  finite  number  of  optimisation 
parameters  (OPs)  to  be  solved  by  usual  meth¬ 
ods.  Because  of  their  convergence  robustness, 
DMs  give  the  safest  approach  for  rapid 
prototyping  of  spatial  trajectories  for  an  a/c  [9]. 

Although  DMs  ensure  much  better  robust¬ 
ness  in  comparison  with  indirect  methods,  it  is 
still  impossible  to  use  them  onboard  for  on-line 
optimisation.  That  is  why  in  expectation  of  more 
powerful  computers,  it  was  proposed  to  employ 
the  databases  of  initial  guesses  of  OPs  or  the 
“banks  of  trajectories”  (TBs)  essentially  im¬ 
proving  a  convergence  robustness. 


Pilot’s  Associate 


“Safety’ 


T 


“Failure” 


“Doctor” 


Information-control  environment 


I  Pilot  J —  - 


A 

A 


E  “Navigator”  J 

1  PSS  -  pilot’s  control  actions  during  maneuvering  support  system 

F  “Maneuver”  ■ 

|]  takeoff 

■ 


takeoff  and  climbing  H  mrface-based  targets  attac  m  |  descent  and  landing 


|  flight  on  a  route 


aircombat 


Fig.3  General  structure  of  PA  expanded  with  PSS 


The  left  and  right  “boundaries”  of  the  tra¬ 
jectory  are  in  general  full  defined.  If  any  of 
states  is  not  fixed,  it  can  be  chosen  with  the  help 
of  corresponding  knowledge-bases,  created  on 
the  basis  of  modelling  and  pilots’  interviewing 
(see  example  on  Fig.5).  Controls’  constraints  are 
supposed  to  be  known  from  the  other  PA’s  sub¬ 
systems. 

I  Choice  of  surface-based  target  attack  direction  I 


The  present  paper  discloses  some  basic  al¬ 
gorithms  of  the  proposed  PSS  and  is  organised 
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as  follows.  Section  2  provides  a  mathematical 
foundation  of  one  of  modifications  of  the  men¬ 
tioned  DM  for  rapid  prototyping  (DMRP)  and 
shows  some  simulation  results.  Section  3  briefly 
discusses  semi-natural  experiments  and  flight 
test  performed  at  that  point. 

2  Real-time  trajectory  optimisation 

2.1  General  statement  of  optimisational 
problem 

The  most  general  statement  of  the  optimal  con¬ 
trol  problem,  determining  a/c  trajectory  from  the 
current  position  to  a  given  point,  may  be  speci¬ 
fied  as  follows. 

There  is  a  set  of  admissible  trajectories 
satisfying:  i)  the  system  of  ordinary  differential 
equations,  ii)  initial  conditions  and  final  (termi¬ 
nal)  conditions,  iii)  constraints  on  the  state 
space,  on  controls,  and  on  the  controls’  deriva¬ 
tives.  The  problem  is  to  find  an  optimal  trajec¬ 
tory  that  minimises  some  cost  function  (CF)  and 
an  optimal  control  corresponding  to  this  trajec¬ 
tory.  Mentioned  CF  in  general  can  be  repre¬ 
sented  not  only  as  integral  function  (the  simplest 
examples  are  the  manoeuvre  time  or  fuel  con¬ 
sumption),  but  also  as  the  function  of  current 
co-ordinates  and  controls  in  the  terminal  point 
or  at  some  event-conditioned  instant  of  time 
(e.g.,  terminal  load  factor  or  bank  angle  at  aim¬ 
ing  point,  etc.). 

Although  the  terminal  point  in  the  preced¬ 
ing  definition  is  considered  as  completely  de¬ 
fined,  in  general,  some  state  variables  and/or 
controls  might  not  be  pre-set.  In  this  case,  they 
are  also  needed  to  be  optimised. 

As  a  system  of  ordinary  differential  equa¬ 
tions,  the  3D  point-mass  equations  over  a  flat 
Earth  with  zero  sideslip  angle  are  considered  as 
usual  (totally  seven  equations).  In  this  case  the 
common  controls  are  a  throttle,  load  factor  and 
bank  angle  with  appropriate  constraints.  Non¬ 
linear  aerodynamic  characteristics  and  propul¬ 
sion  performance  are  presented  by  correspond¬ 
ing  coefficients  defined  by  multidimensional 
tabulated  data. 


2.2  Direct  methods  in  flight  dynamics 

It  was  already  mentioned  that  the  direct  methods 
are  the  best  base  for  the  rapid  trajectory 
prototyping.  The  main  idea  of  the  direct  meth¬ 
ods  is  to  consider  a  function  as  a  finite  set  of 
variables.  This  is  fairly  evident,  if  it  is  assumed 
that  the  admissible  function  can  be  represented 
by  an  infinitive  power  series,  or  by  a  Fourier 
series,  or  by  any  series  of  the  form 


y(x)  ~^lak(Pk(x)’  where  <pk (x)  are  given 

k=\ 


functions.  Then  CF  will  be  the  function  of  a  set 
of  unknown  coefficients  allowing  to  reduce  the 
task  by  considering  a  finite  series  instead  of  in¬ 
finite. 

The  first  direct  approaches  to  solve  the 
problem  of  integral  minimisation  were  intro¬ 
duced  by  Euler  and  later  by  Ritz.  Another  more 
simple  but  more  universal  procedure  was  intro¬ 
duced  by  Galerkin  in  1915  as  a  means  of  ob¬ 
taining  approximate  solutions  to  the  boundary 
value  problems.  When  combined  with  interpo¬ 
lation  equations  of  the  method  of  finite  ele¬ 
ments,  which  is  a  variation  of  the  Rayleigh-Ritz 
procedure,  Galerkin’ s  method  becomes  a  very 
useful  procedure  for  solving  both  initial  and 
boundary  value  problems. 

It  was  Taranenko,  who  developed  and  ap¬ 
plied  the  method  like  Ritz-Galerkin  to  the 
problems  of  flight  dynamics  with  constraints  on 
state  variables  and  controls  [10].  Following  the 
main  idea  of  the  direct  methods,  Taranenko 
suggested  defining  the  reference  functions  (RFs) 
for  both  a/c’s  Cartesian  co-ordinates  and  its  air¬ 
speed  as  x  =  xi0  +  — - —  (t-t0)+<F.(t), 


Tf  T0 


i  =  1,4.  Here  t  is  an  argument,  and  <3>;(t)  are 


continuous,  unequivocal  and  differentiable 
functions,  satisfying  to  obvious  boundary  con¬ 
ditions  dv  (t0  )  =  dx  (tj  )  =  0 .  T aranenko  pro¬ 
posed  to  use  one  of  the  following  functions 

52  At  sin  kn  —  — ,  or  ^Ak{z -tJ(t -rf)  , 

ti  T/-T  o  f~i 

or  (t-t0  Ht-t,)*,  or  their  linear  combina- 
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tion.  However,  in  principle,  there  are  no  limita¬ 
tions,  so  one  can  use  any  convenient  functions 
for  the  particular  task  at  hand. 

An  explicit  mean  to  increase  flexibility  is 
to  increase  a  number  of  elements  n  in  a  men¬ 
tioned  series  or  to  increase  powers  ml  and  m2 . 
However,  Taranenko  proceeded  with  another 
approach  -  he  subdivided  an  interval  [r0;Ty  ]  on 

several  pieces  and  employed  low  order  polyno¬ 
mials  to  describe  behaviour  of  state  variables 
x  ,  i  =  1,4  along  each  of  them  (the  parameters 
of  pieces’  collocation  can  then  be  considered  as 
an  additional  OP).  The  higher  number  n  (or  m] 
and  m2 ),  and  the  higher  the  number  of  pieces  in 
piecewise  case,  the  closer  a  near-optimal  (NO) 
solution  is  to  the  optimal  one. 

The  choice  of  an  argument  T  also  depends 
on  a  particular  task.  Generally  speaking,  one  can 
use  any  continuous  monotonic  parameter:  time, 
path,  etc.  In  case  of  defining  both  a/c’s  co¬ 
ordinates  and  airspeed,  as  Taranenko  did,  we 
may  use  only  an  abstract  parameter,  otherwise 
we  would  be  unable  to  vary  trajectory  and  speed 
history  independently.  Taranenko  called  T  a 
virtual  arc  [10]. 

Finally,  to  be  able  to  satisfy  boundary  con¬ 
ditions  imposed  on  state  variables,  their  first  and 
second  derivatives,  Taranenko  and  his  followers 
preferred  the  sum  of  three  cubic  polynomials  as 
RFs,  where  polynomials’  coefficients  are  auto¬ 
matically  defined  by  pre-set  boundary  condi¬ 
tions  at  the  origin  and  the  end  of  the  trajectory 

[11] .  In  the  case  of  control’s  constraint  viola¬ 
tion,  they  shift  to  a  direct  integration  of  state 
equations  with  the  marginal  value  of  this  (vio¬ 
lated)  control.  Since  the  60’ s,  with  the  use  of 
this  method,  a  lot  of  optimisation  problems  have 
been  solved  for  different  a/c,  including  helicop¬ 
ters,  strike  and  civil  a/c,  airspace  vehicle. 

Neljubov  uses  fifth-seventh  order  time- 
polynomials  for  Cartesian  co-ordinates  only  and 
implemented  them  as  “basic  trajectories”  during 
the  following  tracking  with  a  real  automation 

[12] .  Hargraves  [13],  Convay  [14],  and  others 
use  the  so-called  collocation-based  or  direct 
transcription  methods,  which  is  similar  in  many 


respects  to  Galerkin’s  procedure.  They  reduce 
the  initial  problem  by  segmenting  the  time  in¬ 
terval  into  the  5-20  pieces  and  representing  the 
solution  both  for  state  variables  and  controls  by 
piecewise  polynomials  (constants).  The  tens  of 
unknown  coefficients  are  then  determined  by 
enforcing  continuity  at  the  nodes  and  by  satis¬ 
fying  the  differential  equations  at  some  specified 
points  in  each  segment.  Seywald  [15]  and  others 
use  the  so-called  differential  inclusion  approach. 
They  eliminate  controls  from  the  state  equations 
by  employing  a  description  of  the  dynamical 
system  in  terms  of  its  attainable  set.  Lu  [16] 
uses,  for  each  piece,  an  approach  similar  to  Ta¬ 
ranenko’ s  method.  He  calls  it  inverse  dynamics 
approach.  For  optimisation  of  planar  trajectory 
for  aerospace  vehicle  at  each  of  the  twenty 
pieces  he  pre-determines  one  of  the  state  vari¬ 
ables’  time  history  and  one  of  the  controls’  time 
history  by  cubic  splines,  and  then  he  solves  the 
inverse  task  of  flight  dynamics. 

Unfortunately  neither  of  preceding  meth¬ 
ods  has  properties  required  by  PSS  concept. 
Ones  approaches  do  not  provide  analytical  rep¬ 
resentation  of  the  entire  trajectory,  another  im¬ 
ply  relatively  difficult  numerical  calculations 
with  numerous  OPs.  The  accuracy  of  piecewise- 
based  approaches  directly  depends  on  the  num¬ 
ber  of  segments  used  in  the  approximation.  Fi¬ 
nally,  although  all  cited  approaches  have  been 
successfully  used  for  off-line  trajectory  optimi¬ 
sation,  nobody  tried  to  employ  DM  for  on-line 
optimisation. 

2.3  Introduction  to  the  rapid-prototyping  ap¬ 
proach 

Without  loss  of  generality  let  us  consider  a 
computational  procedure  of  the  basic  modifica¬ 
tion  of  DM,  suitable  for  short-term  spatial  ma¬ 
noeuvres  like  takeoff/climbing,  curved  ap¬ 
proach,  etc.  [17].  This  modification  combines  a 
number  of  advantages  over  both  previously 
mentioned  methods  by  Taranenko  and  Nel¬ 
jubov,  and  consequently  a  close  position  to  Lu’s 
approach.  Though  having  less  possibilities  of 
varying  the  trajectory  itself,  this  algorithm  as¬ 
sures  the  following:  i)  the  boundary  conditions 
are  satisfied  ‘a  priori’;  ii)  an  a/c  control  is  physi- 
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cal  and  realisable  (smooth),  meaning  a  pilot  can 
easily  perform  it;  iii)  the  iterative  process  con¬ 
verges  well,  making  it  possible  to  proceed  with 
on-line  optimisation;  iv)  RT  has  an  analytical 
representation. 

We  take  the  RFs  for  an  a/c’s  co-ordinates 
xi ,  i  =  1,3  as  algebraic  polynomials  of  degree  n 


with  the  virtual  arc  z  as  an  argument: 

^  ^  (max(U-2))!r* 

Xi{z)  =  \^aik - - -  (1) 

t=o  k\ 

The  degree  n  of  these  polynomials  is  de¬ 
termined  by  the  number  of  boundary  conditions 
to  be  met,  so  that  all  coefficients  aik  were  de¬ 
termined  algebraically,  rather  than  varied.  The 
higher  the  maximum  degree  of  time  derivative 
of  an  a/c  co-ordinate  at  initial  and  terminal 
points,  whose  values  (the  derivatives)  are 
known,  the  higher  the  degree  of  the  polynomial. 
The  minimum  degree  of  the  polynomial  is 
n  -  d0  +  df  + 1 ,  which  is  greater  by  one  than  the 

sum  of  the  maximum  orders  of  the  time  deriva¬ 
tive  of  the  a/c  co-ordinates  at  the  initial  and  ter¬ 
minal  points  ( d0  and  df,  respectively)  [12]. 

For  example,  if  we  consider  the  task  with¬ 
out  pre-setting  the  initial  and  terminal  values  of 
second  time  derivatives  of  a/c  co-ordinates’ 
(proportional  to  controls),  that  means 
d0-df  -  1 ,  the  minimum  polynomials’  order  is 

n  -  3 .  Substituting  the  corresponding  values  of 
x0 ,  x'j0 ,  i  =  1,3  for  t  =  0 ,  and  xif ,  xif ,  i  =  1,3 


for  z  =  zf  into  Eq.(l)  (where  zf ,  the  length  of  a 

virtual  arc,  is  considered  as  the  first  OP),  we 
obtain  a  set  of  12  linear  algebraic  equations  for 
12  unknown  coefficients  ajk ,  /'  =  1,3 ,  k  —  0,3 


being  resolved  as: 


ai0  ~  Xi0  ’  ai2  ~ 


2  x'  +  Ax' 

if  ~  '^iO 


+  6 


xif  -  xi0 


'  f 


'/ 


_  *r+*.-0  xif~xi0 

ail  ~Xi0’  ai3  ~  6  2  ^  3 


(2) 


yf  yf 

Of  course,  we  can  compound  the  RFs  as 
superposition  of  several  cubic  polynomials 
[10,11],  and  then  satisfy  the  boundary  condi¬ 
tions  for  controls.  Otherwise,  we  should  employ 


the  higher-order  polynomials.  For  instance, 
fifth-order  polynomials  allow  to  satisfy  the 
boundary  values  for  the  a/c  co-ordinates,  their 
first  and  second  time  derivatives  at  both  ends  of 
the  trajectory  (d0  =  df  =  2).  To  have  a  final  part 

of  the  trajectory  more  smooth,  it  is  better  to  em¬ 
ploy  a  case  when  df  =  3  with  x”  =  0 ,  i  =  1,3 . 

The  only  OP  so  far  was  a  length  of  a  virtual 
arc.  However  sometimes  to  make  RT  more 
flexible,  it  is  worth  to  add  some  more  OPs.  For 
instance,  we  can  add  one  fictive  boundary  con¬ 
dition  x^ ,  i  =  1,3  to  the  case  d0  =  2,  d  f  =  3 , 
and  for  n  =  l  obtain  relations  for  24  coefficients 


aik,  i  =  1,3,  k  =  0,7  in  the  same  manner,  as  in 


(2).  Now  we  can  use  these  fictive  boundary  val¬ 
ues  as  additional  OPs. 

Since  RT  is  defined  not  in  the  time  frame, 
it  does  not  explicitly  determine  a  history  of 
speed.  This  gives  a  great  advantage  -  an  a/c  can 
fly  along  the  same  trajectory  with  different 
speed  histories.  In  general,  the  dependence  V(z) 
may  be  determined  either  by  pre-setting  a  sepa¬ 
rate  RF  V(z),  as  in  Taranenko’s  method,  or  by 


integrating  the  corresponding  state  equation  for 
a  speed  with  predetermined  thrust  history 


V'(z)  =  g(nx  -  sin  y) 


dt  g(nx-sm.Y) 


(3) 


dz  X(z) 
where  A  =  f  is  a  virtual  speed,  7  is  a  flightpath 


angle,  and  nx  is  a  tangential  projection  of  a  load 

factor.  It  means  that  we  can  explicitly  employ 
the  results  of  controls’  synthesis  obtained  with 
the  help  of  indirect  methods.  We  will  further 
deal  specifically  with  the  last  approach,  assum¬ 
ing  that  throttle  versus  time  (arc)  history  is 
known  qualitatively  beforehand. 

For  example,  for  the  time-optimum  prob¬ 
lem,  the  optimum  thrust  control  is  the  on/off 
control.  Hence,  when  solving  such  optimisation 
problem,  we  can  set  several  (for  the  majority  of 
cases  at  hand  -  two)  switching  points:  from 
maximum  thrust  8T  to  minimum  ST  .  at  the 

1  iiid.x  1  min 


moment  z*  ( t*T ),  and  back  from  8T min  to  8T max 


at  the  moment  t”  (  fT*)  (0  <z*T  <z*T*  <zf). 
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Therefore,  the  search  of  the  NO  thrust’s  control 
will  be  made  among  three  admissible  thrust 
histories.  ST  max  -  8r  nini .  ST  niin  -  8r  max  and  8r  max  - 

8t  nilI1  -  8t  max .  Of  course,  for  other  type  of  CF  we 

can  try  other  reasonable  thrust  time-histories. 

In  practice  we  only  have  the  boundary  val¬ 
ues  of  state  variables  and  controls,  so  we  need  to 
calculate  boundary  conditions  for  the  derivatives 
of  the  state  variables  with  respect  to  the  argu¬ 
ment  r  down  to  the  order  of  d0  and  df  re¬ 
spectively  (2).  We  can  make  it  with  the  help  of 
known  kinematic  and  dynamic  equations  of  an 
a/c,  using  obvious  relations  x\  =  X1xi  and 

x"i=k~2[xi-xiX'\,  i  =  1,3  [11].  Corresponding 
values  of  A  and  A'  at  the  boundary  points  are 
determined  as  A0;/  =  ,  X0.f  =  V0.fV^f  . 

During  numerical  solution,  the  parameters 
of  RT  are  calculated  in  N  points  equidistantly 
placed  over  the  virtual  arc,  so  that 
At  =  zf  (N  - 1)”1 .  This  sampling  period  for 

j  - 1  ,N  - 1  corresponds  to  the  time  intervals 

=2^L(t/--T.h)2(V,  +vj-iY  •  With  these 

values  of  At  and  At.,  the  virtual  speed  A  is 
calculated  at  each  step  according  to 
A;.  =  ATAt  T1 . 

The  explicit  laws  for  a/c’s  co-ordinates  (1), 
with  account  of  calculated  speed  history  (3), 
uniquely  determine  the  a/c  attitude  -  flightpath 
and  azimuth  angles,  -  and  remaining  controls  - 
bank  angle  (j)  and  normal  projection  of  load 

factor  nz .  It  can  be  done  through  the  solution  of 
the  inverse  problem  of  flight  dynamics  (have  to 
define  the  controls’  time-histories  corresponding 
to  a  desired  RT). 

Therefore,  for  each  set  of  OPs  we  calculate 
the  value  of  a  CF  along  with  the  value  of  a  pen¬ 
alty  function,  summarising  possible  violations 
of  constraints  and  terminal  speed  error.  As  a  re¬ 
sult,  we  reduce  the  original  problem  to  a  non¬ 
linear  programming  problem,  meaning  we  ob¬ 
tain  a  problem  of  minimisation  for  the  scalar 


function  of  several  (but  not  tens  as  in  Refs.  13- 
lb)  OPs. 

Because  of  erroneous  gradient  information 
(because  of  tabulated  character  of  aerodynamic 
and  thrust  data,  because  of  event-conditioned 
step-changing  mass  or  aerodynamic  configura¬ 
tion,  etc.),  zero-order  algorithms  are  preferred  to 
quadratic  programming.  Another  and  possibly 
the  most  important  reason  for  the  use  of  these 
simple  algorithms  is  that  since  we  are  going  to 
implement  them  on  board  of  an  a/c,  we  need  a 
probability  of  solution  equalling  to  one,  mean¬ 
ing  an  absolute  reliability.  In  practice,  a  number 
of  major-loop  iterations  required  by  any  non¬ 
gradient  algorithm  to  converge  the  task  did  not 
exceed  25-30  iterations  (with  an  arbitrary  initial 
guess).  With  account  of  searching  iterations,  the 
total  number  of  CF  evaluations  was  an  average 
of  100-120.  The  run  time  tcpu  for  IBM486-type 

processor  was  not  more  than  3%  of  trajectory 
duration  itself. 

2.4  Examples  of  near-optimal  manoeuvring 

By  now  within  the  PSS  paradigm  with  the  help 
of  proposed  DMRP,  tens  of  thousand  trajecto¬ 
ries  have  been  calculated  and  tested  for  different 
types  of  a/c  for  such  stages  of  flight  as  take¬ 
off/climb,  surface-based  target  attack  (SBTA) 
and  curved  landing  approach  (CLA).  Figs.6-8 
demonstrate  some  examples  of  such  trajectories, 
calculated  for  the  multi-regime  fighter  with  the 
use  of  seventh-order  polynomials. 

Fig.6  illustrates  the  possibility  of  applying 
DMRP  to  calculate  SBTA  (a)  and  CLA  (b)  tra¬ 
jectories.  Fig.7  shows  an  example  when  the  fi¬ 
nal  states’  manifold  is  determined  by  the  termi¬ 
nal  speed,  terminal  range  to  the  origin  of  the  in¬ 
ertial  frame  and  by  final  flightpath  angle,  but  the 
final  azimuth  angle  is  not  pre-set.  In  this  case, 
the  final  azimuth  angle  was  also  optimised. 

Since  DMRP  takes  into  account  the  pre¬ 
scribed  values  of  high  order  derivatives  at  the 
boundary  points  there  is  no  problem  to  calculate 
(optimise)  a  series  of  manoeuvres  ensuring  a 
smooth  change  of  controls.  Fig.8  gives  an  ex¬ 
ample  of  such  SBTA  trajectory  composed  of 
two  manoeuvres  optimised  one  after  another 


661.7 


Oleg  A.  Yakimenko 


(the  terminal  state  of  the  first  manoeuvre  serves 
as  an  initial  state  for  the  second). 

a) 


V — c — V==5 


clx 


Note,  this  feature  allows  implementation  of 
the  same  approach  for  “dynamic”  optimisation, 
when  the  constraints  on  the  state  variables  or 
desired  final  condition  change  in  time.  The  latter 


means  that  a  trajectory  should  be  recalculated 
from  the  current  (prognosis)  condition  every 
several  seconds,  as  it  can  take  place  during  col¬ 
lision  avoidance  in  free  flight  or  evasion- 
pursuing  in  a  dog-fight  (Fig.9),  for  example. 


Fig.8  Illustration  of  consistent  manoeuvres  calculation 


Fig.  10  illustrates  the  example  of  NO  pur- 
suing-evading  in  an  aggressive  dogfight  with  the 
use  of  DMRP  by  both  sides.  As  a  CF,  the  kill 
probability  computed  with  the  use  of  original 
approach  [18]  was  used  in  this  case.  Another 
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feature  of  a  dogfight  is  that  a  final  position  of  an 
a/c  and  its  velocity  is  not  known  beforehand. 
Therefore  they  were  considered  as  additional 
OPs  (in  total  nine  OPs  were  used  for  each  side). 

As  mentioned  earlier,  DMRP  has  a  very 
important  advantage  from  the  standpoint  of  em¬ 
ploying  it  in  PSS  -  the  rapidity  of  two-point 
boundary-value  problem  solution  even  by  means 
of  non-gradient  method  (tCPU  =  0.03 tf ).  To  cor¬ 
rectly  account  a  run  time,  the  predicted  on  tcpu 

a/c’s  state  should  be  used  as  an  initial  point  for 
RT  calculation.  For  example,  for  Pentiumll- 
class  processors  it  means  no  more  than  one- 
second  prognosis  for  60-s  manoeuvre.  Unfortu¬ 
nately  today’s  onboard  computers  are  much 
slower  than  Pentium-type  ones.  Should  we  wait 
until  the  last  one  will  be  available  on  board,  or  it 
is  possible  to  employ  the  proposed  DMRP  right 
away?  The  answer  is  ‘yes’,  we  can  use  it  on  to¬ 
day’s  a/c,  but  we  have  to  employ  a  known  idea 
of  initial-guess  database. 

2.5.  Trajectories  Databases  Employment 

The  idea  of  convergence  robustness  improve¬ 
ment  through  the  employment  of  databases  of 
initial  guesses  is  fairly  clear.  If  some  kind  of 
TB,  containing  the  values  of  OPs  for  a  certain 
number  of  node  trajectories,  would  be  available 
on  board  computer,  then  we  could  use  it  to  ob¬ 
tain  an  initial  guess  of  OPs  for  the  current  two- 
point  boundary-value  problem  (see  Fig.  11a). 
The  developed  DMRP  allows  creating  these 
TBs  due  to  two  main  reasons.  First,  the  entire 
trajectory  is  defined  only  by  a  few  parameters. 
Second,  these  parameters  have  a  clear  physical 
sense,  and  they  change  smoothly  with  a  change 
of  boundary  conditions  and  constraints’  rigidity. 

Fig. lib  demonstrates  an  example  of  effec¬ 
tiveness  of  TB  employment.  Suppose  we  have 
two  node  trajectories  7j  and  T2 ,  defined  by  sets 

of  OPs  Sj  and  32,  and  differ  from  each  other 
only  by  state  x20.  Then  as  initial  guess  for  an 
arbitrary  trajectory  with  x20  e  [v^ ;  J ,  we  can 
use  the  vector  S*  =  (1- p)E°pt  +  pZ°2pt , 


where  p  — 


.  As  it  turns  out,  a  trajectory 


calculated  with  the  use  of  S*  and  a  trajectory 
finally  optimised  with  the  use  of  this  initial 
guess  practically  coincide,  meaning  that  the  in¬ 
terpolated  vector  S*  is  fairly  close  to  the  opti¬ 
mised  one  Eopt  (this  means  that  S*  is  a  good 
guess  for  3 vpr ). 


0  20  40  60 

Time  (s) 

Fig.ll  Illustration  of  TB  usage 


In  general  case,  there  is  more  then  one  in¬ 
put  parameter;  therefore  a  multi-parametrical 
(multi-entry)  interpolation  is  used.  The  only 
question  is  what  set  of  variables  should  be  used 
as  an  entry  parameters,  and  how  many  nodes  for 
each  of  them  should  be  established.  Obviously, 
the  more  nodes  for  each  entry  parameter  TB  has, 
the  better  the  convergence  robustness  is.  How¬ 
ever  in  general,  there  are  too  many  entry  pa¬ 
rameters:  initial  and  final  values  of  state  vari¬ 
ables,  controls  constraints,  a/c  and  its  engine 
characteristics,  atmospheric  conditions.  As  per¬ 
formed  research  shows,  the  DMRP  is  more  sen¬ 
sitive  to  changes  of  only  13  of  those  parameters 
and  it  is  enough  to  have  from  two  to  six  nodes 
to  improve  the  convergence  of  algorithm  by  the 
factor  of  three  (e.g.,  for  an  IBM486-type  proces¬ 
sor  it  means  tCPU  <1%).  It  was  found  that  the 
values  of  the  initial  and  terminal  azimuth  angles 
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are  the  most  influential.  For  this  reason,  TBs 
should  have  the  largest  mesh  point  frequency  for 
these  parameters,  whereas  for  the  majority  of 
other  parameters  it  was  sufficient  to  have  only 
two  nodes:  for  minimum  and  maximum  possible 
(expected)  values.  Consequently,  the  required 
RAM  volume  to  keep  the  OPs  for  all  varieties  of 
those  entry  parameters  is  reasonable.  For  in¬ 
stance,  the  TB  for  onboard  computation  of 
SBTA-type  trajectories  using  an  IBM486-type 
processor  contains  47,040  trajectories  and  cov¬ 
ers  (with  respect  to  the  target)  any  initial  and 
final  azimuth  angles,  initial  and  final  velocity  in 
the  range  of  [170;250]m/s',  initial  range  up  to 
15,000m,  final  range  within  [l,800;3,600]m, 
initial  altitude  within  [200;l,000]m,  final  diving 
angle  within  [10;30]°,  any  operative  constraints 
on  controls.  To  store  the  values  of  three  OPs  and 
the  value  of  CF,  it  requires  less  then  1Mb  RAM. 
For  other  stages  of  flight  like  CLA  or  take¬ 
off/climb  because  of  a/c  restriction  to  a  runway, 
the  number  of  TB’s  node  trajectories  is  ap¬ 
proximately  eight  times  less. 

Another  good  thing  about  TB  -  it  ensures  a 
unit  probability  of  the  NO  trajectory  computa¬ 
tion  for  a  certain  time  due  to  its  construction. 
Moreover,  before  implementation  on  board,  it 
can  and  has  to  be  cross-checked  in  a  series  of 
intermediate  points. 

2.6  Waypoint  flight 

A  little  bit  different  modification  but  with  the 
use  of  the  same  basic  principles  was  designed 
for  the  solution  of  the  tasks  with  “soft”  (pene¬ 
trative)  restrictions  on  state  co-ordinates  like  a 
waypoint  flight  through  the  field  of  prohibited 
zones  of  plight  (PZF)  (air  threat  zones,  thunder¬ 
storm  fronts,  etc.)  [19]. 

In  this  modification  instead  of  pre-setting  a 
thrust  history,  the  separate  RF  for  speed  was 
employed.  Following  the  basic  idea  of  DM,  the 
search  of  NO  trajectories  was  made  in  the  class 
of  global  cubic  splines  passing  through  the  set 
of  waypoints,  location  of  which  have  to  be  op¬ 
timised.  Fig.  12  explains  the  essence  of  this 
modification. 


Fig.12  Illustration  of  spline-variations  method  essence 

Generalised  function  of  losses,  compound 
of  accumulated  damage  and  a  total  fuel  con¬ 
sumption  during  entire  flight  (with  coefficient 
k),  was  considered  as  a  CF  for  this  kind  of  ma¬ 
noeuvres.  The  real  values  of  coefficient 
k  s  [l0“5;10“3]  factually  lead  to  neglecting  of 
the  functional  losses,  but  in  case  of  lack  of  PZF 
we  automatically  receive  a  NO  solution  for  the 
fuel-minimum  flight.  Because  of  polymodal 
nature  of  the  CF  in  this  task,  Strongin’s  infor¬ 
mation-statistical  method  [20]  was  employed  in 
order  to  find  an  area  of  global  extremum  attrac¬ 
tion.  To  assure  the  required  robustness  a  zero- 
order  method  was  employed  then  for  a  final  lo¬ 
cal  optimisation. 

Fig.  13  shows  an  example  of  a  route  opti¬ 
misation  in  presence  of  PZF  field,  formed  by 
five  stationary  air  threat  zones,  at  variation  of 
different  number  of  states  (free  and  fixed  speed 
history).  Terrain  features  are  also  accounted 
here.  Fixing  one  of  the  states  (some  concrete 
tasks  might  require  certain  altitude  or  speed  pro¬ 
file),  we  can  receive  different  sub-optimal  solu¬ 
tions.  It  demonstrates  a  flexibility  of  the  de¬ 
signed  DM. 

Fig.  14  demonstrates  one  more  example  of 
spatial  trajectory  optimisation  in  condition  of 
flight  through  the  complex  PZF  field,  which  be¬ 
sides  three  fixed  circle-type  zones  and  one  tri- 
angle-type  zone,  contains  also  a  moveable  PZF, 
crossing  the  line  of  pre-set  a/c  path  with  a  con¬ 
stant  speed. 

Optimisation  algorithm  turned  to  be  fairly 
sensitive  to  a/c  performance  (it  smartly  uses  the 
features  of  a  concrete  a/c  dynamics).  Where  a 
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condition  of  non-growth  of  damage  permits,  a/c 
flies  on  fuel-optimum  regime. 
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Fig.13  Examples  of  NO  waypoint  flight 


Fig.  14  Route  optimisation  in  presence  of 
moveable  PZF 


flight  simulator  (SpMATS)  on  the  modern 
multi-regime  fighter  [21]. 

SpMATS  represents  a  complex  consisting 
from  the  real  a/c,  local  network  of  two  field  PC 
linked  to  onboard  computer,  and  FD  (Fig.  15). 
First  PC  carries  out  a  modelling  of  an  a/c’s  6- 
DOF  dynamics  (an  object-oriented  approach 
allows  to  quickly  change  an  a/c  type),  synthesis 
of  cockpit’s  instrumental  panel,  HUD  images 
and  synthetic  outcabin  world,  as  well  as  PSS 
prototype  functioning.  SpMATS  functioning 
requires  an  a/c  to  be  under  current  with 
switched-on  objective  control  system  (OCS)  and 
released  control  column  (without  start  of  en¬ 
gines).  The  pilot,  sitting  inside  the  a/c  cabin, 
with  the  help  of  FD  observes  synthetic  world, 
targets,  basic  indicators,  and  HUD  images  cor¬ 
responding  to  a  current  flight  mode  and  simula¬ 
tion/training  task.  He  operates  with  the  a/c  and 
armament  control  system  models  by  usual  con¬ 
trols.  Modelling  PC  takes  these  controls  inputs 
from  onboard  PC  (OCS). 


PC  of  modeling  PC  of  control 

Fig.15  SpMATS  block-scheme 


Run  time  for  this  modification  depends  on 
task  dimensionality,  making  in  average  of 
tcpu  =  5-10%  for  IB M466-type  computer. 

3  Testing  and  onboard  implementation  expe¬ 
rience 

Proposed  PSS’s  algorithms  were  tested  not  only 
in  PC  simulation,  but  also  in  real  simulation 
with  pilots. 

3.1  Research  simulator 

A  semi-natural  experiment  was  mainly  done 
with  the  use  of  adoptive  cockpit-transformed 


An  instructor/researcher  uses  the  second 
PC  for  the  management  of  simulation/training 
process.  He  can  not  only  supervise  the  training 
process,  but  also  participate  in  it  as  instructor 
(co-pilot)  or  opponent.  He  also  can  introduce 
any  disturbances,  a/c’  systems  failures  and  Ac¬ 
tive  damages,  new  targets,  etc.  Both  PC  are 
linked  through  the  standard  bi-directional  port 
RS-232  and  use  a  specially  developed  protocol. 

LCD-based  FD  is  installed  with  the  help  of 
rubber  bandage  in  front  of  HUD  (Fig.  16). 
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3.2  Semi-natural  experiments 

During  testing  with  the  help  of  SpMATS  pilots 
were  asked  to  observe  different  NO  RT,  which 
were  generated  on-line  for  real  tactical  condi¬ 
tions  (Fig.  17  shows  examples  of  FD  formats). 
RiS  was  equipped  with  appropriate  wayside 
signs,  like  “Throttle-up”,  “Throttle-down”, 
“Fire”,  “Drag  flaps”,  “Gears”,  etc. 


Fig.16  FD  accommodation 


Fig.17  FD  formats:  a)  SBTA,  b)  waypoint  flight 


In  total  14  pilots  participated  in  the  semi¬ 
natural  experiment  on  PSS’  algorithms  testing 
[22],  Some  estimates  and  recommendations 
were  also  given  by  test-pilots  Pugatchev, 
Bichkov,  Kutuzov,  Votintzev,  Solovjov,  and 
Tzoj.  Totally  there  were  analysed  over  a  hun¬ 


dred  and  fifty  realisations  of  different  manoeu¬ 
vres.  For  each  type  of  manoeuvres  an  accuracy 
of  trajectory  observation,  an  accuracy  of  termi¬ 
nal  conditions  satisfaction,  a  character  of  pilot’s 
control  actions,  and  dependence  on  flight  regime 
were  analysed.  The  PSS’s  and  RiS’s  parameters, 
e.g.  RiS’s  cross-section’s  size,  RiS’  type  (gutter, 
tunnel);  image  type  (real,  compressed),  etc.  were 
analysed  and  optimised  also.  Of  course  all  ex¬ 
periments  were  carried  our  with  the  model  of 
real  atmosphere. 

The  vigorous  (high-G)  trajectories,  specifi¬ 
cally  SBTA  and  pursuing-evading  trajectories, 
as  expected,  were  the  most  difficult  ones  for 
tracking.  However  even  for  them,  the  accuracy 
of  RT  tracking  was  good  enough:  I D  ~  60 m , 

/AV  ~-10 mis,  /|An  |  =0.6,  /|Ar|  -10°,  AT  ~  3s 
(here  h  =  \  f  ('))*< dt,  T  is  a  du- 

1  T  V  i= 2 

ration  of  manoeuvre,  IM  =  —  j  (a(t)  -  aref  (t))dt ). 

T  T 

Terminal  conditions  were  also  satisfied  with  a 
high  accuracy  (final  attitude  errors  in  vertical 
and  horizontal  planes  were  of  order  of  two  de¬ 
grees). 

Experiments  proved  that  RTs  can  be  effec¬ 
tively  calculated  on-line,  and  that  RiS-image 
keeps  all  necessary  information  for  their  effec¬ 
tive  tracking.  It  was  also  confirmed  that  the  pi¬ 
lots  master  a  new  control  mode  with  the  RiS 
assistance  fairly  quickly.  Fig.  18  shows  one  ex¬ 
ample  when  after  the  proper  training  the  pilot 
could  perform  a  manoeuvre,  which  is  closed  to 
the  NO  one,  himself,  without  prompt  (it  means 
PSS  paradigm  may  be  used  also  for  pilots’ 
training).  The  spectral  analysis  of  control  ac¬ 
tions  proves  pilots’  subjective  opinion  that  there 
is  no  additional  workload  because  of  this  new 
kind  of  indication. 

Obviously,  the  accuracy  of  tracking  de¬ 
pends  upon  RiS’s  cross-section  dimensions 
(Fig.  19)  and  RT  curvature.  That  is  why  for  the 
CLA  trajectories  the  averaged  tracking  accuracy 
was  much  better  -  I D  ~  10m . 
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The  final  errors  at  glide-path  capture  made 
the  order  of  5 m/s  in  terminal  speed,  ~  4m  in 
lateral  deviation,  ~2m  in  altitude,  ~2°  and  ~1°  in 


terminal  pathangle  and  azimuth.  Obviously  they 
did  not  affect  on  the  subsequent  gliding  and 
landing. 
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Fig.18  SBTA  trajectory  tracking 
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Fig.19  Accuracy  and  intensity  of  pilot's  control  actions 


3.3  Flight  test 

Finally,  designed  software  was  flight-tested  for 
the  CLA  trajectories  by  the  test  pilots  Birjukov 
and  Kazin  onboard  the  flying  laboratory  An 
(Antonov)-72  equipped  with  Glonass  and  Nav- 
star/GPS  data  receivers  in  the  Gromov  Flight- 
Test  Institute  [23],  FD  was  placed  in  front  of 
pilots  over  the  instrument  panel.  On-line  RT 


computation  was  provided  onboard  by  the  sepa¬ 
rate  IBM486-type  computer. 

It  was  shown  that  proposed  FD  with  PSS’ 
algorithms  allowed  performing  CLA  with  the 
high  accuracy  and  pilot’s  confidence.  Pointing  a 
prognosis  marker  along  the  RiS  provided  the 
movement  inside  the  corridor  by  light  adjusting 
control  actions.  Operating  by  the  throttle,  the 
pilot  provided  the  required  decrease  of  airspeed. 
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Several  CLA  were  performed  in  twilight  condi¬ 
tion  and  with  strong  wind  turbulence.  Fig.20 


shows  an  example  of  final  part  of  CLA  in  con¬ 
dition  of  lateral  wind  transition  of  10mA. 


Fig.20  Example  of  final  straight-in  approach  with  lateral  wind 


Concluding  Remarks 

In  conclusion  it  can  be  stated  that  the  developed 
DM  solve  optimisational  boundary-value  prob¬ 
lems  very  effectively.  Resulting  algorithms  that 
have  an  excellent  convergence  robustness  can 
be  easily  integrated  with  existing  naviga¬ 
tion/control  algorithms,  making  possible  the  on¬ 
line  prototyping  of  spatial  manoeuvres  for  their 
following  tracking  by  a  pilot  with  the  use  of  FD. 
The  high  efficiency  and  necessity  for  pilots  of 
such  kind  of  support  was  convincingly  proved  in 
ah  carried  out  semi-natural  experiments  and 
flight  test.  Therefore,  performed  research  gives 
the  grounds  for  further  proceeding  with  imple¬ 
mentation  of  PS  S  as  a  part  of  PA  onboard  an  a/c 
in  order  to  support  the  lowest  levels  of  pilot’s 
decision  making. 
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